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Abstract 

We demonstrate that a measurement of the impact parameter in one-prong r 
decay can be useful for the determination of the Higgs boson parity in the H/A —* 
t + t - ; — > p~^v T (v T } decay chain. We have estimated that for a detection set-up 
such as TESLA, use of the information from the r impact parameter can improve 
the significance of the measurement of the parity of the Standard Model 120 GeV 
Higgs boson to ~ 4.5a, and in general by factor of about 1.5 with respect to the 
method where this information is not used. 

We also show that the variation in the assumption on the precision of the mea- 
surement of the impact parameter and/or 7r's momenta does not affect the sensitivity 
of the method. This is because the method remains limited by the type of twofold 
ambiguity in reconstructing the r momentum. 
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1 Introduction 



One of the most important goals in the scientific program of a future linear electron- 
positron collider (LC) is to measure the properties of the Higgs boson precisely. Among 
them, the parity plays a prominent role. Depending on the mass and the Higgs mechanism, 
the Higgs parity can be measured by different means [HEllSllllE] • In some cases, it can be 
best determined by investigating the properties of the Higgs boson decay into r leptons [HI 
[7j. In jH] we studied the possibility of distinguishing a scalar from a pseudoscalar coupling 
of a light boson to fermions using its decay to a pair of r leptons and their subsequent 
decays r ± — ► p ± j? T (u T ) and p ± — > 7r ± 7r . This decay chain gives an independent test of 
model and production mechanism. In fact, we have found that the angular distributions of 
the decay products, which clearly distinguish the different parity states, are measurable 
using typical properties of a future detector at an e + e~ linear collider. A special Monte 
Carlo technique was developed for that purpose 0. In the present paper we will extend 
that analysis, using the information that can be obtained from the measurement of the r 
lepton impact parameter. 

It turns out that the kinematic configuration of e + e~ — > t + t~ events cannot be fully 
determined since the momenta of both v T and v T are unknown. The lost information can 
be recovered partially by measuring energy and directions of all visible particles [TTH ITT]. 
When both r leptons decay hadronically and all hadron momenta are recovered, the 
original r direction can be obtained up to a twofold ambiguity [1011121 ■ In the method, 
the requirement that the two r leptons be essentially of equal energy and back to back 
in some known frame, is decisive. The measurement of the tracks of the hadrons, in 
particular their relative impact parameters, allows us to resolve this twofold ambiguity, 
see for instance jTUHHllISl • 

In the case of linear colliders such as TESLA [H|, a sufficiently precise reconstruction 
of the frame, where r leptons are expected to be back to back turns out to be impossible, 
mainly due to beamstrahlung effect. This is a strong limitation for the Higgs boson parity 
measurement as proposed in Refs. [37j. A way around was found in Ref. Here, we will 
investigate if the measurement of the r lepton impact parameter direction, with precision 
as expected for TESLA, can be used to optimize that method. 

This work is organized as follows: in section 2 the geometry of r production and decay 
is described, as well as our way to use the r impact parameter in the measurement of the 
Higgs boson parity. In section 3 some details of the TAUOLA Monte Carlo simulation are 
given. Section 4 is devoted to the definitions used in our simulation of detector properties. 
In section 5 definitions of our observable and numerical results are given. The summary, 
section 6, closes the paper. 
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2 Description of the method 



In the decay of the r lepton into one charged particle, the r impact parameter can be 
defined as the distance of closest approach (in the plane perpendicular to the beam axis) 
of the charged particle to a reference point, which is assumed to be the production point 
of the r (see Fig. QJ. A positive sign is given to this quantity if the track crosses the 
r direction in the same hemisphere, relative to the r production point, in which it lays, 
and negative otherwise. The production point is assumed to lie on the beam line. In our 
study we will exploit the possibility to determine the direction of the impact parameter, 
rather than the related distance. In fact we will be interested, where the r production 
point is localized with respect to the charged pion track. 

Before we address the question of the information available from the impact parameter, 
let us recall some details of the actual observable proposed in Ref. [H] for the Higgs boson 
parity measurement. The method relies on measuring the acoplanarity angle of the two 
planes, spanned on decay products and defined in the p + p~ pair rest frame; some 
additional selection cuts also need to be applied. Let us present the elements of the 
observable as the following points. 




Figure 1: Schematic view of the r decay to p and v T . PP is the r 'production point, DP 
the r decay point, and PC A the point of closest approach. Directions and energies ofir, 
tt° , p^ are directly measurable, the energy of the r ± can be reconstructed from the event 
topology. The r momentum must be placed on the cone (drawn on the plot) and the r 
impact parameter can be used to determine the actual position. 



1. Acoplanarity of the p+ and p decay planes: 
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Figure 2: The circle of candidates for r momentum. Depending on the possible geometrical 
configuration, three cases are possible: A) tt direction inside the circle: single candidate 
for r momentum. B) tt direction outside the circle: two candidates for r momentum. 
C) tt direction outside the circle and no crossing of dashed line with the circle (due to 
detection ambiguities of measured angles and energies). We take in this case the closest 
point to the line. Effectively, we obtain a single solution. 



• The reconstructed four-momenta of tt + and tt° are combined to yield the p + 
four- momentum. The same is done for the p~ . 

• All reconstructed four-momenta are boosted into the p + p~ rest frame. 

• The angle (p between the planes of the p + and p~ decay products in this frame 
is the acoplanarity. 

2. Normalized energy differences: 

The events are divided into two classes depending on the value of yiy2, where 

^ _ E w + — E^o _ ^ _ E n - — E^o 
E n + + E n o E v - + E n o 

The energies of tt^, tt° are to be taken in the respective rest frames. 

3. Replacement r rest frame: 

In Ref. the method of reconstruction of replacement rest frames was proposed. 
We will use this method here as well, but we will omit the details. In fact they are 
rather unimportant. The energies of pions taken directly from the laboratory frame 
can equally well be used in the above definitions of Eq. (JTJ). 

4. Higgs rest frame: 

If the information on the beam energies and the energies of all other observed 
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particles is taken into consideration, the Higgs rest frame can be reconstructed, 
for instance in the Higgsstrahlung production process, e + e" — > HZ, when the Z 
boson decays either into a charged lepton pair or hadronically. We may define 
the 'reconstructed' Higgs boson momentum as the difference of the sum of beam 
momenta and momenta of all visible particles, i.e. decay products of Z and all 
radiative photons of | cos#| < 0.98. 

5. r energy: 

• In the reconstructed Higgs boson rest frame, the r four-momenta are estimated 
in a crude way by assuming the direction of the respective p's and an energy 
of mjj/2. 

• The t momenta are boosted back to the laboratory frame and their energies 
are taken to be the reconstructed energies of r leptons. 

6. r direction: 

The r direction in the laboratory system is constrained by two requirements: 

• It has to lie on a cone around the p direction with opening angle ip defined as: 

2E T E p -ml -ml 
cos^j = — S 2 

ZPrPp 

where E, p, m are the energy, momentum and mass, respectively. The four- 
momentum of p is measured as the sum of the four-momenta of its charged 
and neutral daughters. 

• It has to lie in the plane spanned by the vector pointing from PP to PC A and 
the 7r + momentum. 

The intersection of the cone and the plane is calculated numerically. There are three 
cases: 

• One solution: the solution is taken as the r direction (see Fig. |2] case A). 

• Two solutions: one of the two solutions (see Fig. Incase B) is taken on a random 
basis. 

• No solution: the direction on the cone closest to the PCA-n + momentum 
plane is taken (see Fig. El case C). 
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7. Impact parameter improved replacement T-rest frame: 

With the help of the r energy and r momentum defined in points 5 and 6 above, 
a new impact parameter improved replacement r rest frame can be defined. In this 
way we have an alternative way of estimating the difference of 7r ± , 7r° energies see 
Eq. (0) in rest frames. 

3 Monte Carlo set-up 

In the following discussion all the Monte Carlo samples have been generated with the 
TAUOLA library [T5lll6pi7| . For the production of r lepton pairs, the Monte Carlo program 
PYTHIA 6 . 1 [18 was used 1 . The production process e + e" — > ZH — > n + ji~{qq)H has been 
chosen, with a Higgs boson mass of 120 GeV and a centre-of-mass energy of 350 GeV. 

The effects of initial-state bremsstrahlung were included in the PYTHIA generation. 
For the lepton pair decay (with full spin effects from the H — > t + t~, r ± — > p ± u T (u T ), 
p ± — > 7T ± 7r° chain) the interface explained in Ref. was used. It is an extended version 
of the interface of Refs. dUE]]. 

4 Detection parameters 

To test the feasibility of the measurement, some assumptions about the detector effects 
have to be made. We include, as the most critical for our discussion, effects due to inaccu- 
racies in the measurements of the 7r ± , 7r° momenta and r ± leptons impact parameters. We 
assume Gaussian spreads of the 'measured' quantities with respect to the generated ones, 
and we use the following algorithm to reconstruct the energies of 7r's in their respective 
rest frames. 

1. Charged-pion momentum: We assume a 0.1% spread on its energy and direction. 

2. Neutral-pion momentum: We assume an energy spread of , 5% For the 

1 toJ v y/E[GeV] 

9 and (ft angular spread we assume |^q. (These neutral pion resolutions can be 
achieved with a 15% energy error and a 27r/1800 direction error in the gammas 
resulting from the tt° decays.) These resolutions have been approximately verified 
with a SIMDET [22J, a parametric Monte Carlo program for a TESLA detector [23J. 

1 It was shown that the interface can work in the same manner with the HERWIG fj] generator. 
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Figure 3: The acoplanarity distribution angle <p of the p + p~ decay products in the rest 
frame of the p + p~ pair. Detector smearing is included. Generator level r ± rest frames are 
used. The thick line corresponds to a scalar Higgs boson, the thin line to a pseudoscalar 
one. The left figure contains events with yiy 2 > 0, the right is for y x y 2 < 0. 

3. The reconstructed Higgs boson rest frame: We assume a spread of 2 GeV 
with respect to the transverse momentum of the Higgs boson, and 5 GeV (because 
of the beamstrahlung effect) for the longitudinal component. 

4. The impact parameter: 

The angular resolution of the vector pointing from PP to PC A (see Fig. [1} has 
been simulated for the TESLA-like detector. The simulation is based on the an- 
ticipated performance of a 5-layer CCD vertex detector as described in For 
Higgsstrahlung events with H —>■ t + t~ and r ± — ► p ± u T ( y u T ) at run = 120 GeV and 
■\/s = 350 GeV, the angular resolution has been found to be approximately 25°. We 
use this resolution in our simulations. 
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Figure 4: The acoplanarity distribution angle <p of the p + p~ decay products in the rest 
frame of the p + p~ pair. Full smearing is included. Replacement r^ 1 rest frames defined as 
in Ref. J3j/ are used. The thick line corresponds to a scalar Higgs boson, the thin line to a 
pseudoscalar one. The left figure contains events with y^y 2 > 0, the right is for y x y 2 < 0. 

5 Acoplanarity of the p + and p~ decay products with 
the help of the impact parameter 

For all figures in our paper we will use the same angle <p built from smeared four-momenta 
of the p ± decay products. The usefulness of this acoplanarity distribution manifests itself 
only after selection cuts are applied. 

In Fig. El our p + p~ decay products acoplanarity distribution angle <p, defined in the 
rest frame of the reconstructed p + p~ pair is shown. Unobservable generator-level rest 
frames are used for the calculation of selection cuts (Eq. ((TJ)). The two plots represent 
events selected by the differences of rr^ tt° energies defined in these respective r ± rest 
frames. On the left plot, the signs are required to be the same yiy 2 > 0, whereas on the 
right one, events are taken with yiy 2 < 0. This figure quantifies the size of the parity effect 
in an idealized condition, which we will attempt to approach with realistic ones. The size 
of the effect was substantially diminished when the detector set-up was included for 
rest frames reconstruction as well, see Fig. 01 following exactly the method of Ref. [Sj. 
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Figure 5: The p + p~ decay products' acoplanarity distribution angle, (p, in the rest frame 
of the p + p~ pair. Full smearing is included. The r lepton impact parameter is used in the 
reconstruction of the t^ 1 rest frames. The thick line corresponds to a scalar Higgs boson, 
the thin line to a pseudoscalar one. The left figure contains events with yiy 2 > 0, the right 
is fory x y 2 < 0. 

Let us use this result later, as a reference point to quantify the size of our improvements. 

It turns out that a slight improvement is visible, when the impact parameter is taken 
for the rest frames determination, as explained in section 4. In fact, the maximum 
relative difference between the scalar and pseudoscalar cases is larger by ~ 12% (compare 
Figs. EJand EI) if a Gaussian spread of the r impact parameter o~ip = 25° is used. We 
have also checked that the precision of the r lepton impact parameter measurement is not 
critical in this case. The twofold ambiguities as explained in Fig. 121 are the main reason 
for the sensitivity to be much lower than in Fig. El 

At the cost of introducing cuts and thus reducing the number of accepted events, we 
can go one step further in the improvement of the method. We first start by requiring 
that the signs of the energy differences y\ and y 2 be the same with both methods (with 
and without r lepton impact parameter). The result is depicted in Fig. El Only ~ 52% 
of events are accepted, but the parity effect grows by ~ 107% in comparison to Fig. |U 

The events can also be treated differently, according to the number of solutions; see 
Fig. El If the same cut is imposed only for leptons with double solutions, while for 
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Figure 6: The p + p~ decay products' acoplanarity distribution angle, (p, in the rest frame 
of the p + p~ pair. Only events where the signs of the energy differences yl and y2 are the 
same whether calculated using the method described in Ref. JSJ/ or with the help of the r 
impact parameter are taken. The thick line corresponds to a scalar Higgs boson, the thin 
line to a pseudoscalar one. 

single-solution events the r rest frames obtained with the help of the impact parameter 
are used alone, the number of accepted events turns out to be ~ 72%, but the parity 
effect is larger by only ~ 67% than in Fig. HJ This result is depicted in Fig. 

Other cuts are also envisageable, but the two proposed here are easy to implement and 
lead to an acceptable loss of statistics, together with a sizeable increase of parity effect so 
that there is not much room for further improvements. 



6 Summary 

We have found that the measurement of the impact parameter can be helpful in the 
determination of the Higgs boson parity at a future linear collider, using the H/A — > t + t~; 

— > p ± z/ r (z/ T ) cascade decay. Thanks to the improved method the effect of parity is 
enhanced by ~ 107% at the cost of losing only ~ 48% events, which represents a gain 
of a factor of ~ 1.5 in sensitivity with respect to the analysis without impact parameter. 
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Figure 7: /ie p + p~ decay products' acoplanarity distribution angle, (p, in the rest frame of 
the p + p~ pair, with cuts depending on the position of the pion with respect to the r cone. 
The selection is as in Fig. G| but in the case of single- solution events (Fig. [H case A), the 
result from the old method is ignored. The thick line corresponds to a scalar Higgs boson, 
the thin line to a pseudoscalar one. 

This correspond to ~ 4.5cr separation between scalar and pseudoscalar for the 120 GeV 
Standard Model Higgs and other assumptions similar to the one taken in Ref. [S|. 

We have also checked, that reasonable variations in the assumptions of detection un- 
certainties in measurement of charged and neutral 7r's as well as of the r + and r~ impact 
parameters, hardly change the sensitivity of our method at all. This is because the dom- 
inant uncertainty is saturated by the twofold ambiguity as explained in Fig. 2. We are 
afraid that a reconstruction of the Higgs boson momentum to precision better than the r 
mass would be needed to resolve the ambiguity. This seems unrealistic for several reasons, 
e.g. the beamstrahlung effect. 

Probably, a fit to the multidimensional distribution might be more efficient than the 
simple cuts applied in this work. We believe, however, that such studies should be done 
only when the detector properties have been understood better. Finally, we expect that 
further improvements, e.g. involving the r ± — > afu T decay mode, may turn out to be 
helpful for the measurement of the Higgs boson parity. 
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